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Abstract: The plantation of non-native species is one of the most expensive ecological restoration measures 
in arid and semi-arid areas, while its impacts on local communities are largely ignored. This study assessed the 
rate of change and the dynamic degree of the economic values of ecosystem services related to local 
conservation (water yield, stocking rate and aesthetic value) and preserving the future (carbon sequestration, 
soil protection, soil stability and habitat provision) to determine the restoration success of the plantation of 
non-native species Haloxylon ammodendron (C.A.Mey.) Bunge ex Fenzl (15- and 30-year-old) in parts of arid 
rangelands of Bardsir region, Kerman Province, Iran. We investigated the impacts of the two plantations on 
the seven ecosystem services and ecosystem structures (horizontal and vertical structures, vegetation 
composition and species diversity) based on field sampling and measurements at four sampling sites (i.e., 
control, degraded, and 15- and 30-year-old plantation sites) in spring and summer of 2022. The restoration 
success of the plantation of non-native species was then examined by assessing the rate of change and the 
dynamic degree of the total economic value of all ecosystem services as well as the rate of change and the 
dynamic degree of the economic values of ecosystem services for the two groups (local conservation and 
preserving the future). Although the plantation of non-native species H. ammodendron enormously improved 
the vertical and horizontal structures of ecosystems, it failed to increase species diversity and richness fully. 
Further, despite the plantation of non-native species H. ammodendron had significantly increased the economic 
values of all ecosystem services, it was only quite successful in restoring carbon sequestration. Path analysis 
showed that plantation age had a significant impact on restoration success directly and indirectly (through 
changing ecosystem structures and services). The dynamic degree of the economic values of ecosystem 
services related to local conservation and preserving the future at the 15- and 30-year-old plantation sites 
indicated that the two plantations successfully restored the ecosystem services related to preserving the 
future. The presented method can help managers select the best restoration practices and predict their 
ecological-social success, especially for the plantation of high-risk non-native species in arid and semi-arid 
areas. 
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1 Introduction 


Arid and semi-arid areas cover approximately 42.00% of the Earth's total land area, supporting 
38.00% of the world's population, who are often the poor ones (Huang et al., 2017). Land 
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degradation in arid and semi-arid areas is accelerating, and ecosystems are facing serious 
exploitation problems (Delgado and Marin, 2020). Inappropriate management plans not only 
cannot lead to improvement, but also may exacerbate land degradation in arid ecosystems 
(AbdelRahman, 2023). The utter ecological methods should be reconsidered by decision-makers 
to provide a comprehensive approach to ecosystem management (Ceccon et al., 2020). 
Sustainable land management is necessary for the health of both ecological and social systems 
(Lébmann et al., 2022). However, people and the economic benefits they derive from the 
ecosystems are often ignored in the decision-making process (Chee, 2004). New integrated 
approaches are needed to bridge the gap between ecological and social assessments of ecological 
restoration (Löfqvist et al., 2022). Restoration is one of the ultimate goals of ecosystem 
management, which aims to return degraded ecosystems to the pre-degradation state (Saco et al., 
2006). Restoration may change ecosystem structures (Zhou et al., 2020) and functions (Santos et 
al., 2021). The paucity of information on the impact of these changes is an obstacle to ecosystem 
restoration (Liniger et al., 2011). Moreover, the revision of approaches in monitoring plant 
communities' structures and functions is essential for the future decisions of land management 
(Emamian et al., 2021). 

The plantation of plant species adapted to harsh environments is one of the approaches to 
recovering degraded ecosystems (Vlasenko et al., 2022). Assessment of restoration success is 
very important for land management due to the high cost of restoration projects (Yıldız et al., 
2022). The attributes related to biodiversity (Farahat and Linderholm, 2012), ecosystem structures 
(Yang et al., 2022) and ecosystem functions (Beldini et al., 2010) are measured to assess the 
success or failure of plantation projects. New approaches based on ecosystem services are 
developed for ecological-social assessment of ecosystems in recent years (Peng et al., 2023; 
Sharafatmandrad and Khosravi Mashizi, 2023). Sustainable restoration projects should provide 
benefits for both ecological and social systems by improving ecosystem services (Bullock et al., 
2011; Keesstra et al., 2018). The economic value of ecosystem services reflects the value of 
nature (Millennium Ecosystem Assessment, 2005), showing not only the economic importance of 
ecosystems (Sannigrahi et al., 2020) but also their health status (Yan et al., 2016). Assessing the 
economic value of ecosystem services under land cover change is an appropriate tool, which can 
help decision-makers and politicians protect ecosystems and contribute to sustainable 
development (Liu et al., 2018b). The relationship between land cover change and ecosystem 
services has been studied widely (Reyers et al., 2009; Ouyang et al., 2016; Gao et al., 2017). 
Plantation may decrease or increase the supply of ecosystem services by changing vegetation 
composition (Reyers et al., 2009; Ouyang et al., 2016; Zhang et al., 201 6b). 

Plant species that conserve ecological integrity and support the well-being of social system 
should be considered for plantation in arid and semi-arid areas (Reisman-Berman et al., 2019). 
Both native and non-native species are commonly used for plantation in these areas. However, the 
plantation of native species does not change ecosystem structures (Wang et al., 2022), while the 
plantation of non-native species may strongly alter ecosystem structures, functions and processes 
due to their different growth forms compared to the native ones (Vu Ho et al., 2023). Previous 
studies have shown both positive and negative impacts of planting non-native species on 
ecosystems (Bravo et al., 2019; Randriambanona et al., 2019). The plantation of non-native 
species can alter the biodiversity of ecosystems by changing soil properties, litter decomposition 
rate, microclimate conditions, etc. (Baker and Murray, 2012). Planting non-native species with 
thick trunks may reduce the abundance of native species (Harrington and Ewel, 1997). However, 
non-native species may have positive effects on soil function (Zhu et al., 2020). Randriambanona 
et al. (2019) showed that afforestation by pine as a non-native species had acted as a catalyst and 
increased the abundance of native species in the ecosystem. Unfortunately, most studies only 
focused on tree species and their edaphic impacts (Hoque et al., 2021), timber production and 
carbon sequestration (Ye Myint et al., 2021), while ignoring other ecosystem services. In addition, 
monoculture restoration planting is rapidly growing across the world (Liu et al., 2018a). 
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Reduction in species richness is one of the challenges of monoculture with non-native species 
(e.g., Rédei et al., 2020; Vu Ho et al., 2023). Planting non-native species requires a 
comprehensive ecological assessment in order to understand the possible impacts 
(Reisman-Berman et al., 2019). There is a paucity of information on the impacts of planting 
non-native species on ecosystem services in arid and semi-arid areas. The plantation of non-native 
species for ecological restoration needs to proceed with caution, especially in fragile arid 
ecosystems. 

Therefore, a comprehensive approach is needed to evaluate the success of restoration projects 
in order to assist managers in decision making. In addition to ecological benefits, restoration 
projects should provide sufficient economic benefits to locals (Zeng et al., 2021). Monitoring the 
changes of plant communities and the economic value of ecosystem services provides a 
theoretical reference and a basic scientific method for assessing the success of restoration 
projects. Current methods of assessing the success of restoration projects often tend to focus on 
the ecological aspects and neglect the social aspects. In this study, we assessed the variations in 
seven ecosystem services related to local conservation (water yield, stocking rate and aesthetic 
value) and preserving the future (carbon sequestration, soil protection, soil stability and habitat 
provision), as well as ecosystem structures (horizontal and vertical structures, vegetation 
composition and species diversity), to determine the restoration success of 15- and 30-year-old 
plantations of non-native species Haloxylon ammodendron (C.A.Mey.) Bunge ex Fenzl in parts of 
arid rangelands of Bardsir region, Kerman Province, Iran. The study pursues the following 
objectives: (1) the effect of non-native species plantations on ecosystem structures; (2) the effect 
of non-native species plantations on the economic values of ecosystem services; (3) the 
relationship between ecosystem structures and ecosystem services; and (4) the assessment of 
restoration success through assessing the rate of change and the dynamic degree of the economic 
values of ecosystem services related to local conservation and preserving the future at the 15- and 
30-year-old plantation sites. The comprehensive method presented in this study is a step ahead of 
the previous methods to assess the restoration success of the plantation of non-native species from 
an ecological-social point of view. 


2 Materials and methods 


2.1 Study area 


This study was carried out in parts of arid rangelands of Bardsir region (29°40'-30°05'N, 
56°40'—-57°20'E; 1862-2300 m a.s.l.), Kerman Province, southeastern Iran. The landscape is an 
alluvial plain with deep sandy loam soils that are loose and highly susceptible to wind erosion. The 
region covers a total area of 9.34x10* hm’, with a mean annual precipitation of 256 mm (which 
often occurs between October and December). The long-term mean annual temperature is 15.6°C, 
the mean maximum temperature is 35.0°C (which usually occurs in June and July) and the mean 
minimum temperature reaches —10.0°C (which usually occurs in January and December). 

The dominant plant species of this region is Artemisia sieberi Besser., which is severely 
degraded due to human overexploitation. Non-native species H. ammodendron had been planted 
to restore the ecosystems in the region. H. ammodendron is a Mediterranean species that is 
widespread in Asia and Africa. Haloxylon genus is highly adapted to arid and semi-arid areas and 
can grow well in areas with annual precipitation of less than 100 mm (Wickens et al., 1985). 
Therefore, this genus is widely used to restore arid and semi-arid areas and stabilize sand dunes 
due to its high resistance to drought, salinity and heat (Arokh et al., 2021). H. ammodendron can 
obtain water sources differing in space and time by its dimorphic root system (Dai et al., 2022). 
Branches of H. ammodendron have photosynthesis capability, indicating the adaptation of this 
species to the low humidity and high temperature in drylands (Hu et al., 2021). In Iran, about 
2.00x10° hm? of land area has been afforested by Haloxylon sp., especially in central and 
southeastern regions of the country (Arokh et al., 2021). 
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2.2 Sampling and measurements of ecosystem structures 


In the study area, four sampling sites were selected: control, degraded, and 15- and 30-year-old H. 
ammodendron plantation sites. Sampling sites were similar in terms of topography, soil types and 
geological formations. At the 15- and 30-year-old H. ammodendron plantation sites, H. 
ammodendron plantation projects have been implemented in 2007 and 1992, respectively, to 
restore vegetation and reduce wind erosion. At the control site, grazing intensity is low (Holechek 
and Galt, 2000) and the late successional plant species A. sieberi is dominant. At the degraded 
site, the cover of palatable and late successional plant species is reduced due to overgrazing and 
land use change, and thus annual species Salsola brachiata Pall. is dominated. 

At each sampling site, the horizontal and vertical structures of ecosystems, vegetation 
composition and species diversity were measured as ecosystem structures. Landscape Function 
Analysis (LFA) method was used to evaluate the vertical and horizontal structures of ecosystems. 
This method was applied in areas with different climate types from dry rangelands in Australia 
(Tongway and Smith, 1989) to rainforests in Indonesia (Tongway and Hindley, 2003). The 
accuracy of this method for conditions in Iran has been confirmed by various studies (e.g., Ata 
Rezaei et al., 2006; Dehghani Bidgoli and Keshavarz, 2018; Sharafatmandrad and Khosravi 
Mashizi, 2019). For this method, landscapes are divided into patches (perennials and litter) and 
inter-patches (bare open soil) (Tongway and Hindley, 2004). Three randomly selected 100 mx100 
m plots were used to collect data at each sampling site. A 50-m transect was laid out in the 
direction of the prevailing wind in each plot. The relative area (%) of patches and inter-patches 
was measured along each transect. Landscape Organization Index (LOI) and canopy volume 
(m°/hm?) of patches were considered as the horizontal and vertical structures of ecosystems, 
respectively. LOI reflects the proportion of an ecosystem that is able to capture and utilize vital 
resources (Tongway and Hindley, 2003); it can be calculated by the proportion of the length of 
patches to the total length of transect (Tongway and Hindley, 2003). Mean canopy area 
(m*/hm7) and height (m) of patches were used to measure the canopy volume of patches in 
different height classes (Tongway and Hindley, 2003; Saaed et al., 2022). Statistical analyses of 
data on the length, width and height of patches and the canopy volume in different height 
classes (Table 1) were performed using LFA method (Tongway and Hindley, 2003). All field 
sampling and measurements were conducted on May and June of 2022. 


Table 1 Height classes of patches (perennials and litter) based on Landscape Function Analysis (LFA) method 


Height class Height (m) Height class Height (m) 
1 0.00-0.50 4 1.50-2.00 
2 0.50-1.00 5 2.00-2.50 
3 1.00-1.50 6 2.50-3.00 


Simpson's diversity index was used to reflect species diversity as follows: 
D= Soy (1) 
i=l 


where D is the Simpson's diversity index; s is the plant species; and p; is the proportion of the 
number of i" species divided by the total number of all plant species. 


2.3 Quantification of the value of ecosystem services 


In this study, seven ecosystem services related to local conservation and preserving the future 
were quantified (Table 2). The Integrated Valuation of Ecosystem Service and Tradeoff (InVEST) 
model was applied to quantify water yield, carbon sequestration and soil protection. Soil stability, 
aesthetic value and habitat provision were calculated using field data. Metabolizable energy in 
forage and metabolic energy required by livestock were applied to estimate the stocking rate at 
each sampling site (Coupland, 1992; Robles and Passerat, 1995). The detailed description on 
quantifying the seven ecosystem services is provided in Appendix. 
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The market price method was used to measure the economic values of carbon sequestration 
(Ninan and Inoue, 2013), water yield (de Groot et al., 2002) and stocking rate (Baldwin et al., 
2022), since there were market prices for these services. The contingent valuation method (CVM) 
was applied to assess the economic values of non-market services, i.e., soil stability, aesthetic 
value and habitat provision. CVM is widely used for estimating non-market values (de Groot et 
al., 2002) and is a stated preference (survey) method that directly asks people about their 
willingness to pay for a particular ecosystem service (Gürlük, 2006). Previous studies showed that 
CVM is useful in measuring the economic values of soil stability, habitat provision (de Groot et 
al., 2002) and aesthetic value (Garcia-Llorente et al., 2012; Rewitzer et al., 2017) of ecosystems. 
It should be noted that in this study, we converted the economic value of Iranian Rials (IRR) into 
the international currency (USD) based on the exchange rate in 2022 (1 USD=42,000 Rials). 


Table 2 Ecosystem services quantified in this study 


Service category Service type Reference 
ot ak : Stocking rate 
Provisioning services f Havstad et al. (2007) 
Water yield 


; 7 Carbon sequestration 
Regulating services . . Pan et al. (2013) 
Soil protection 
Cultural service Aesthetic value van Zanten et al. (2016) 
Habitat provision de Groot et al. (2002); Millennium 


Supporting services Soil stability Ecosystem Assessment (2005) 


2.4 Determination of restoration success 


In this study, we divided the selected seven ecosystem services (water yield, stocking rate, carbon 
sequestration, soil protection, soil stability, habitat provision and aesthetic value) into two groups 
based on management goals: local conservation and preserving the future (Himes et al., 2020). 
Water yield, stocking rate and aesthetic value are services that local people are willing to protect, 
because these services are directly related to the well-being of locals. Regulating and supporting 
services that are important for the sustainable production of services in the ecosystem must be 
preserved as intermediary services for future ecosystem sustainability (Millennium Ecosystem 
Assessment, 2005). Therefore, the ecosystem services of water yield, stocking rate and aesthetic 
value were included in local conservation group and the ecosystem services of carbon 
sequestration, soil protection, soil stability and habitat provision were included in preserving the 
future group. The restoration success of plantations was then examined by assessing the rate of 
change and the dynamic degree of the total economic value of all ecosystem services as well as 
the rate of change and the dynamic degree of the economic values of ecosystem services related 
to the two groups (local conservation and preserving the future). 

The rate of change, which indicates the size and magnitude of the change, was used to measure 
the changes in the economic values of ecosystem services and ecosystem structures at the 
degraded, and 15- and 30-year-old H. ammodendron plantation sites compared to the control site. 
It can be calculated as follows (Liu et al., 2018b): 

AV 


A 
AV =Vr—Va ’ (3) 
where R (%) is the rate of change; AV (USD) is the change in the economic values of ecosystem 
services (or ecosystem structures); and Vr (USD) and V, (USD) are the economic values of 
ecosystem services (or ecosystem structures) at the control site and other sites, respectively. 
The dynamic degree of the economic values of ecosystem services is a quantitative value used 
to compare the economic values of ecosystem services between plantations with different ages 
and predict their future trends (Liu et al., 2018b). It is calculated as follows: 


R= x100% , (2) 
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cate <4 <100% , (4) 
ae 
where K (%) is the dynamic degree of the economic values of ecosystem services; V, (USD) and 
Va (USD) are the economic values of ecosystem services at the plantation and degraded sites, 
respectively; and T (a) is the plantation age. 

Ecosystem structures (horizontal and vertical structures, vegetation composition and species 
diversity) and services (provisioning services, regulating services, supporting services and 
cultural service) as well as plantation age were considered as the drivers of restoration success. 
Path analysis was used to examine the complex relationship between restoration success drivers. 
It is an advanced statistical method by which, in addition to direct effects, the indirect effects of 
independent variables on dependent variables can be identified (Lande and Arnold, 1983). Path 
coefficients can be calculated using the standardized regression coefficients. Regression 
coefficients are obtained by establishing structural equations, i.e., equations that determine the 
structure of the assumed relationships in a model (Fan et al., 201 6b). 


2.5 Statistical analysis 


We used one-way analysis of variance (ANOVA) followed by the least significant difference 
(LSD) to compare the four sampling sites based on data on the economic values of the seven 
ecosystem services (soil stability, aesthetic value, carbon sequestration, water yield, soil 
protection, habitat provision and stocking rate) and ecosystem structures (horizontal and vertical 
structures, vegetation composition and species diversity). Before conducting one-way ANOVA, 
we checked data for the normality and homogeneity of variance. Pearson's correlation was used to 
determine the relationship between ecosystem structures and ecosystem services. 


3 Results 


3.1 Variations in ecosystem structures and ecosystem services 


Canopy volume was distributed in three height classes at the control site, with the maximum 
value of 1200 m*/hm? (Fig. 1). At the degraded site, canopy volume was distributed in three 
height classes, with the maximum value of 320 m?/hm’. At the 15-year-old plantation site, canopy 
volume was distributed in four height classes, with the maximum value of 932 m?/hm?. Further, at 
the 30-year-old plantation site, canopy volume was distributed in six height classes, with the 
maximum value of 928 m3/hm?. 

The relative area of patches with perennials decreased from 41.21% (+1.23%) at the control 
site to 14.32% (+2.31%) at the degraded site (Table 3). The relative area of patches with 
perennials at the 15- and 30-year-old plantation sites were 31.20% (+1.32%) and 44.21% 
(43.12%), respectively. The LOI values of the control and 15- and 30-year-old plantation sites 
were 0.86, 0.68 and 0.73, respectively (Table 3). This means that 86.00% of the area at the control 
site includes patches where resources are conserved, and the other 14.00% of the area at the 
control site consists of inter-patches where resources may be lost from the ecosystem. For the 15- 
and 30-year-old plantation sites, 32.00% and 27.00% of the areas consist of inter-patches where 
resources may be lost from the ecosystem, respectively. 

Subshrubs showed the highest percentage of vegetation cover at the control (27.30% (+3.65%)) 
and degraded (10.25% (41.32%)) sites (Table 4), while small trees exhibited the highest 
percentage of vegetation cover at the 15-year-old plantation (25.23% (+5.36%)) and 30-year-old 
plantation (37.23% (+5.37%)) sites. The highest value of Simpson's diversity index occurred at 
the control site (0.65 (+0.13)). The Simpson's diversity index values were 0.52 (£0.13) and 0.58 
(+0.14) at the 15- and 30-year-old plantation sites, respectively. 

Figure 2 shows the changes in ecosystem services (water yield, stocking rate, carbon 
sequestration, soil protection, habitat provision, soil stability and aesthetic value) at the control, 
degraded, and 15- and 30-year-old plantation sites. Ecosystem services were improved at the 
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30-year-old plantation site compared to the 15-year-old plantation site. With the exception of 
carbon sequestration, all other ecosystem services exhibited the highest and lowest values at the 
control and degraded sites, respectively. Specifically, carbon sequestration showed the highest 


value at the 30-year-old plantation site. 


The economic value of stocking rate at the control site was 329.20 (£42.54) USD/hm?, which 
was higher than the economic values of other ecosystem services (with the exception of carbon 
sequestration; Table 5). The economic values of carbon sequestration were highest at all sampling 
sites among all ecosystem services. With the exception of carbon sequestration, the highest 
economic value of ecosystem services was related to the control site. Moreover, there was no 
significant difference of economic values between the control and 30-year-old plantation sites in 


terms of soil stability (P>0.05). 
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Fig. 1 Distribution of canopy volume in different height classes at the control site (a), degraded site (b), 


15-year-old plantation site (c), and 30-year-old plantation site (d) 


Table 3 Characteristics of horizontal structure at the control, degraded, and 15- and 30-year-old plantation sites 


Horizontal structure Control site Degraded site 


15-year-old plantation site 30-year-old plantation site 


Relative area of Perennials 41.2141.23° 14.3242.31° 31.20+1.32° 44.21+3.12° 
patches (%) Litter 8.32+3.45° 1.34+1.65° 6.23+1.34° 8.56+2.14° 
Relative area of Bare open soil 50.4743.12" 84.3443.12° 62.5743.21° 47.23+3.12° 
inter-patches (%) 
LOI 0.86 0.34 0.68 0.73 


Note: Mean+SD. Different lowercase letters within the same row indicate significant differences among the different sampling sites 


(P<0.05), while the same letter indicates no significant difference. 


Table 4 Vegetation cover and species diversity at the control, degraded, and 15- and 30-year-old plantation sites 


Index Control site Degraded site 15-year-old plantation site 30-year-old plantation site 
Shrubs 3.56+2.40* 1.05+0.56° 1.12+0.56* 2.14+1.03" 
Vegetation Subshrubs 27.3043.65° — 10,25+1.32° 1.23+0.35° 3.13+1.02* 
cover (%) Herbs 4.2340.56"  3.23+0.35° 3.2141.20° 3.24+1.32° 
Small trees 0.00+0.00° 0.00+0.00° 25.23+5.36” 37.23+5.37° 
Species D 0.65+0.13° 0.21+0.06° 0.52+0.13° 0.58+0.14° 
diversity Richness 26+3° 845" 1244? 13+5° 


Note: D, Simpson's diversity index. Mean+SD. Different lowercase letters within the same row indicate significant differences among 


the different sampling sites (P<0.05), while the same letter indicates no significant difference. 


Mohsen SHARAFATMANDRAD et al.: Evaluation of restoration success in arid... 1297 


E 3.0 F (a) 80r (b) 
2 a 
3 x 
% > 
2 15+ = 40 
- S 
2 3 
z| | 
3 2 
n 
0 Control Degraded 15-year-old 30-year-old ° Control Degraded 15-year-old 30-year-old 
plantation plantation plantation plantation 
Site S Site 
T g 
«i 20 Ff (c) < 
E 2 507 @) 
oD = 
Žž 5 
& 10 E 25 
i) & 
2 3 
g a 
a a 
= g 
A Control Degraded 15-year-old 30-year-old E Control Degraded 15-year-old 30-year-old 
plantation plantation 5 plantation plantation 
Site Site 
0.8 f (e) 
2 gr © 
‘3 £ 
2 04 3 15 
o 
v 3 
< 0.0 L L L J > L L L J 
“Control Degraded 15-year-old 30-year-old Control Degraded 15-year-old 30-year-old 
plantation plantation plantation plantation 
0.8 - (2) Site Site 


Habitat provision 
> 
aS 
T 


Control Degraded 15-year-old 30-year-old 
plantation plantation 


Site 
Fig. 2 Variations in ecosystem services of stocking rate (a), soil stability (b), soil protection (c), carbon 


sequestration (d), aesthetic value (e), water yield (f) and habitat provision (g) at the control, degraded, and 15- and 
30-year-old plantation sites 


Table 5 Economic values of ecosystem services at the control, degraded, and 15- and 30-year-old plantation sites 


Economic value (USD/hm’) 


Ecosystem service 


Control site Degraded site 15-year-old plantation site 30-year-old plantation site 

Stocking rate 329.20442.544 18.40+2.32 89.80£13.54° 128.40418.95° 

Water yield 134.03441.57° 49.3848.23* 98.10+14.23° 103.38+9.16° 

Carbon sequestration 653.30+85.21° 228.30+71.70* 445.10+70.50° 827.80+98.204 
Soil protection 12.48+1.41° 6.74+0.75° 8.38+1.05° 11.67+1.06° 
Habitat provision 10.2141.35° 3.21+0.23° 5.35+0.34° 8.35+0.38° 
Soil stability 14.23+0.35° 5.68+0.35* 9.3841.23° 13.2541.54° 
Aesthetic value 18.35+2.68° 4.21+41.23* 10.2543.25° 12.32+4.39° 


Note: Mean+SD. Different lowercase letters within the same row indicate significant differences among the different sampling sites at 
P<0.05 level, while the same letter indicates no significant difference. 
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3.2 Determination of the restoration success of 15- and 30-year-old H. ammodendron 
plantations 


Figure 3 shows the rates of change in horizontal and vertical structures and species diversity at the 
degraded, and 15- and 30-year-old plantation sites compared to the control site. The rates of 
change in horizontal and vertical structures and species diversity were negative at the15-year-old 
plantation site; while they were positive at the 30-year-old plantation site except for the rate of 
change in species diversity. 


40 F $ Vertical structure 4 Species diversity 
20L F Horizontal structure . 
= z 
Ş 0 f fi j 
T I 
v 
2 -20H i 
5 
s —40F 
2 
g ~60F I 
—80F Į 
=100 = 
Degraded 15-year-old 30-year-old 
plantation plantation 


Site 
Fig. 3 Rates of change in horizontal and vertical structures and species diversity at the degraded, and 15- and 
30-year-old plantation sites compared to the control site. Error bars represent standard deviations. 


Comparison of the degraded and plantation sites with the control site showed that the 
maximum negative rate of change was related to the stocking rate at the degraded site (Table 6). 
The rates of change in the economic values of the seven ecosystem services were negative but 
smaller at the two plantation sites than at the degraded site, except for the rate of change of the 
economic value of carbon sequestration at the 30-year-old plantation site, which was positive and 
26.00% higher than that at the control site. In terms of ecosystem services related to the two 
groups (local conservation and preserving the future) and all ecosystem services, the rates of 
change were respectively —48.00% (+3.65%), -36.00% (42.31%), and —42.00% (43.21%) at the 
15-year-old plantation site, and —38.00% (+4.32%), -2.00% (+0.12%) and —20.00% (+5.23%) at 
the 30-year-old plantation site (Fig. 4). 

The dynamic degrees of the economic values of ecosystem services related to the two groups 
(local conservation and preserving the future) and of the total economic value of all ecosystem 
services at the 15-year-old plantation site were higher than those at the 30-year-old plantation site 


(Fig. 5). 


Table 6 Rates of change in the economic values of the seven ecosystem services at the degraded, and 15- and 
30-year-old plantation sites compared to the control site 


Rate of change (%) 


Ecosystem service 


Degraded site 15-year-old plantation site 30-year-old plantation site 

Livestock rate —94.0042.32 —73.004£12.32 —61.00+8.65 
Water yield —63.17+48.21 —27.0048.23 —23.00+£5.68 
Carbon sequestration —65.78+4.23 —31.0045.32 26.00+4.35 
Soil protection —46.0045.65 —32.00+7.65 —06.00+3.20 
Habitat provision —68.00+4.35 —47.00+2.35 —18.00+2.35 
Soil stability —60.00+2.31 —34.00+3.68 —7.00+1.32 

Aesthetic value —77.00+6.23 —44.00+6.98 —32.00+03.68 


Note: Mean+SD. 
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Fig. 4 Rates of change in the economic values of ecosystem services related to the two groups (local 
conservation and preserving the future) and in the total economic value of all ecosystem services at the degraded, 
and 15- and 30-year-old plantation sites compared to the control site. Error bars represent standard deviations. 
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Fig. 5 Dynamic degrees of the economic values of ecosystem services related to the two groups (local 
conservation and preserving the future) and of the total economic value of all ecosystem services at the 15- and 
30-year-old plantation sites compared to the degraded site. Error bars represent standard deviations. 


3.3 Impacts of ecosystem services and ecosystem structures as well as plantation age on 
restoration success 


3.3.1 Correlations between ecosystem structures and ecosystem services 


Pearson's correlation analysis showed that horizontal structure was significantly positively 
correlated with water yield, soil protection and habitat provision (P<0.05). Vertical structure was 
strongly and negatively correlated with stocking rate and positively correlated with carbon 
sequestration and soil protection (P<0.05). Vegetation composition had a strong negative 
correlation with carbon sequestration and a strong positive correlation with stocking rate and 
aesthetic value (P<0.05). Species diversity was positively correlated with stocking rate, water 
yield, habitat provision and soil stability, and negatively correlated with soil protection (Fig. 6). 
3.3.2 Direct and indirect impacts of restoration success drivers 


Using path analysis, we investigated the direct and indirect impacts of restoration success drivers 
(Table 7). The results showed that ecosystem structures, ecosystem services and plantation age all 
had significant impacts on restoration success. Plantation age had the most direct impact on 
restoration success (P<0.05). It also indirectly affected restoration success through ecosystem 
structures and services. The indirect impact of plantation age on restoration success through 
ecosystem structures was greater than that through ecosystem services (Table 8). 
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Horizontal Vertical Vegetation Species 
structure structure composition diversity 


Stocking rate 


Water yield Moderate positive 


Carbon sequestration Strong positive 
Soil protection Non-significant 
Habitat provision Moderate negative 
Soil stability Strong negative 


Aesthetic value 
Fig. 6 Correlations between ecosystem structures and ecosystem services 


Table 7 Direct and indirect impacts of restoration success drivers based on path analysis 


Independent variable Dependent variable Standardized beta t Adjusted R? F 
Model 1 0.85 12.35" 
Plantation age 0.46 13.24" 
Vertical structure 0.34 8.35" 
Horizontal structure 0.25 3.21* 
Vegetation composition —0.37 -9.31** 
Species diversity Restoration success 0.33 7.21" 
Provisioning services 0.23 3.20" 
Regulating services 0.39 10.23" 
Supporting services 0.34 8.35" 
Cultural services 0.27 4.32" 
Model 2 0.83 10.09** 
Vertical structure 0.43 9.35" 
Horizontal structure . 0.30 5.28™ 
j Bes Plantation age 
Vegetation composition 0.39 -7.23 
Species diversity 0.32 6.08" 
Model 3 0.79 7.13" 
Provisioning services 0.37 4.85% 
Regulating services . 0.47 7.28™ 
; ; Plantation age oe 
Supporting services 0.40 5.32 
Cultural service 0.30 3.68" 


Note: ** indicates significance at the P<0.01 level, and * indicates significance at the P<0.05 level. 


Table 8 Standardized impacts of plantation age on restoration success 


Restoration success driver Impact Standardized beta t 
Direct impact 0.46 13.24" 
: Indirect impact through ecosystem structures 0.27 7.67" 
Plantation age i . . n 
Indirect impact through ecosystem services 0.21 6.34 
Total impact 0.90 18.34** 


Note: ** indicates significance at the P<0.01 level. 


4 Discussion 


Non-native small tree H. ammodendron is very different from native subshrub A. sieberi. The 
plantation of H. ammodendron had drastically changed ecosystem structures in arid and semi-arid 
ecosystems in the region. Changes in ecosystem structures will result in variations of ecosystem 
services (Assis et al., 2023). Relationships and processes are as important as structural attributes 
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in ecosystem sustainability (Wu et al., 2022). Thus, the dynamic degree of the economic values of 
ecosystem services can provide very important information for sustainable restoration projects 
and ecological development. 


4.1 Impact of the plantation of non-native species H. ammodendron on ecosystem 
structures 


At the two plantation sites, the plantation of H. ammodendron increased the relative area of 
patches and the height classes of vegetation. There were three more height classes in H. 
ammodendron plantation sites than the control site. Increase of canopy volume in different height 
classes had a very effective role in promoting the vertical and horizontal structures of ecosystems. 
H. ammodendron plants usually grow in large clumps in arid ecosystems and can enhance the 
vegetation cover (Ninot et al., 2007). The larger the clumps, the more reduction the environmental 
disturbances in arid and semi-arid areas (Fryrear, 1995; Marquart et al., 2019). Wind erosion is 
one of the main disturbances in arid and semi-arid areas (Du et al., 2022). H. ammodendron as a 
natural barrier plays an important role in preventing wind erosion through improving the vertical 
and horizontal structures of ecosystems. There are complex plant-soil and plant-plant interactions 
in the ecosystems (Pan et al., 2021). Therefore, the changes in vegetation composition and 
ecosystem functions resulted from the plantation of non-native species should be considered in 
ecosystem management. 

H. ammodendron plantations had led to the formation of large patches in the ecosystems, 
thereby affecting the density and diversity of native species in the study area. Only 37.54% of 
vegetation cover at the degraded site belonged to the dominant subshrubs at the control site. The 
vegetation cover of subshrubs had decreased to 4.50%-—11.46% at the two plantation sites, 
indicating a strong negative impact of H. ammodendron on the dominant native species. Previous 
studies had also shown that the introduction of non-native species into local ecosystems will 
affect the composition of native vegetation (e.g., Kotzen, 2003). Some studies stated that 
monoculture can provide better conditions for native species and lead to higher species diversity 
(e.g., Cuong et al., 2013). Brown (2003) and El-Wahab et al. (2014) indicated that Haloxylon sp. 
could improve plant diversity of predominantly understory annual species under H. ammodendron 
overstory cover. Elmefregy and El-Sheikh (2020) pointed out that soil properties, hydrological 
characteristics, dispersal functions and topographic characteristics are the main environmental 
factors altering the impact of Haloxylon sp. on plant diversity. In contrast, some monocultures 
may prevent the establishment of native seedlings by changing soil properties (Kelty, 2006). H. 
ammodendron affects the composition of native plant communities through changing some 
physical and chemical properties of soil (Jafari et al., 2004). Haloxylon spp. can increase soil 
salinity (due to salt accumulation in the litter layer) and limit the growth of native plant species 
(Rathore et al., 2015). On the other hands, H. ammodendron has higher competitive ability due to 
the development of horizontal and vertical roots that usually inhibit the growth of surrounding 
plants because of allopathic effects (Keneshloo et al., 2018). The higher water absorption capacity 
of H. ammodendron root system is also one of the obstacles to the growth of other plants in arid 
and semi-arid areas. Zhou et al. (2022) reported the gradual reduction of surface soil moisture 
over increasing age of H. ammodendron plantations in China, and concluded that faster soil 
drying threatens the growth of smaller native understory plants. 

At the 15-year-old plantation site in the study area, richness also decreased by 53.85% 
compared to the control site. Natural ecosystems have more diverse vegetation composition and 
age structure than monocultures (Navarro-Cerrillo et al., 2020; Santini et al., 2020). The more 
plant species there are, the greater the diversity of niches; consequently, the accompanying 
understory plant species will be more diverse (Larjavaara, 2008). Plantations are more vulnerable 
to environmental disturbances due to their uniform structures (Camarero et al., 2021). 


4.2 Relationship between ecosystem structures and ecosystem services 


In the study area, species diversity was the most important attribute of ecosystem structures 
influencing the supply of some ecosystem services (stocking rate, water yield, habitat provision 
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and soil stability). Some studies also showed significant positive relationships between species 
diversity and ecosystem services (Quijas et al., 2010; Dee et al., 2017; Himes et al., 2020). Reale 
et al. (2022) pointed out that biodiversity loss strongly affects the global economy through 
changing the supply of ecosystem services. The link between species diversity and ecosystem 
services showed that biodiversity conservation is a desirable economic solution (Dee et al., 2017). 

Provisioning services (water yield and stocking rate) that are directly related to human 
livelihood are more sensitive to the loss of species diversity (Sharafatmandrad and Khosravi 
Mashizi, 2021). Diaz et al. (2006) and Schmid et al. (2009) also indicated that plant diversity can 
increase primary productivity and thus provision of food. Our results showed that supporting 
services (soil stability and habitat provision) that can provide the basis for supplying other 
services are strongly related to species diversity. Higher biodiversity is required for the 
sustainable supply of ecosystem services (Hooper et al., 2005). Increase of species diversity can 
not only increase the diversity of plant use (Pretzsch and Schiitze, 2016), but also contribute to the 
plant communities that are more resilient against environmental stresses such as drought due to 
the complementary effects between species with diverse structures (Gazol and Camarero, 2016; 
Grossiord, 2020). 

Diverse ecosystems provide a wider range of services when moderate level of ecosystem 
functions is desired, but monocultures can maximize the provision of any specific service (van 
der Plas et al., 2016). In the study area, H. ammodendron plantations as monocultures provided a 
relatively long-time protection against wind erosion through improving the vertical structure of 
ecosystems, but they failed in increasing species diversity. This may be caused by the failure of H. 
ammodendron plantations in promoting ecosystem multi-functionality. There are both trade-off 
and synergy relationships between monocultures and the provision of some ecosystem services 
(Suarez and Gwozdz, 2023). Therefore, restoration programs should be aimed at preserving and 
enhancing species diversity to sustain the provision of multiple ecosystem services. Some studies 
showed that plantations with two, three or four plant species are more productive and successful 
in increasing biodiversity (Liu et al., 2018a; Rivera-Pedroza et al., 2019). Although monoculture 
is known as the dominant type of plantations in practice (Rivera-Pedroza et al., 2019), it is not 
different from other diverse ecosystems in providing some ecosystem services. Polley et al. 
(2020) indicated that the stability of biomass production is similar in monoculture and 
mixed-species grasslands. Maintaining natural ecosystems alongside monocultures is considered 
as a strategy to preserve biodiversity in ecosystems where monocultures should be applied 
(Rivera-Pedroza et al., 2019; Yahya et al., 2022). 


4.3 Impact of the plantation of non-native species H. ammodendron on ecosystem services 


The results of this study showed that 15- and 30-year-old plantations respectively enhanced the 
total economic value of all ecosystem services 1.23-5.35 and 2.17—7.76 times, compared to the 
degraded site. Avtar et al. (2014) concluded that plantations can lead to the restoration of 
ecosystem services such as soil protection, carbon sequestration and water conservation, and 
may reduce poverty in developing countries through increasing economic benefits. In the study 
area, the greatest success of plantations was enhanced carbon sequestration (1.8—3.6 times at the 
plantations sites higher than those at the degraded site). Ahmadi kareh et al. (2013), Loni et al. 
(2018) and Ma et al. (2021) also indicated that the plantation of Haloxylon sp. can enhance the 
carbon sequestration potential of arid ecosystems. The results of Kumar et al. (2020) also 
showed that plantations may improve the carbon sequestration of ecosystems in arid and 
semi-arid areas. 

In the study area, ecosystem services were improved at the plantation sites, but this is not the 
only option for sustainable management. What matters in ecosystem management is how much of 
the potential of the control site (as an optimal situation in providing ecosystem services) is met by 
plantations after 15-30 years. Compared to the control site, plantations were only successful in 
providing carbon sequestration service (Table 6). Previous studies also concluded that Haloxylon 
communities can significantly store more carbon than rangelands by enhancing aboveground and 
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belowground biomass (Kareh et al., 2013; Ma et al., 2021). Increasing the carbon sequestration 
potential of ecosystems can modulate the impact of climate change on arid and semi-arid 
ecosystems (Lal et al., 2011), which is very effective in protecting the future of ecosystems and 
mitigating the impact of climate change. 

At the plantation sites, soil protection service was almost slightly less than that at the control 
site (Table 6). Previous studies have shown that soil properties may be improved after planting 
(Singh et al., 2020; Kumi et al., 2021). Soil protection may enhance with increasing silt and clay 
fractions of the soil in areas planted with Haloxylon sp. (Zhang et al., 2016a). Abdi et al. (2019) 
indicated that Haloxylon sp. can increase soil cohesion in arid areas through its root system. Both 
drought and heavy rainfall may threaten plant species in arid and semi-arid areas (Peñuelas and 
Filella, 2001), thereby reducing plant growth rate and biomass in ecosystems (Zamin et al., 2013). 
Soil retention stabilizes plant growth and enhances ecosystem resilience to these threats in arid 
and semi-arid areas (Wolka, 2014; Zhang et al., 2014). In general, plantations were most 
successful in providing regulating services and less successful in providing provisioning and 
cultural services in the study area. Dai et al. (2018) also showed that plantations are more capable 
of recovering regulating services than other services. 


4.4 Impact of plantation age on ecosystem services 


The findings of this study showed that the vertical and horizontal structures of ecosystems are 
improved over increasing plantation age (Fig. 6). Ebrahimi et al. (2017) also indicated that as the 
age of H. ammodendron individuals increases, the cover of plant community is improved, and the 
density and length of the branches increase, thereby forming larger patches. The number of 
services provided by ecosystems may also change over increasing plantation age (Zeng et al., 
2021). 

Regulating, provisioning, supporting and cultural services respectively increased a lot over 
increasing plantation age. As the age of plantation increases, the amount of carbon stored in 
ecosystems increases correspondingly (Zhang et al., 2019; Tamang et al., 2021). Older plants can 
store more carbon due to higher biomass (Fonseca et al., 2012; Wang et al., 2013). In this study, 
plantation age had a greater impact on provisioning services than supporting services. The 
dynamic degree of the economic values of ecosystem services did not increase over increasing 
plantation age. In other words, the dynamic degree of the economic values of ecosystem services 
at the 15-year-old plantation site was higher than that at the 30-year-old plantation site. The rate 
of change in ecosystem services was greater in the younger plantations than in the older 
plantations (Yamaura et al., 2021). 


4.5 Restoration success 


Restoration planting projects are successful when conservation is done by local people. Public 
participation in management programs is usually very limited (Shrestha and Shrestha, 2017; Liu 
et al., 2019). In terms of providing supporting and regulating services that can stabilize 
ecosystems and are necessary for the future of conservation (Millennium Ecosystem Assessment, 
2005), plantations were much more successful in enhancing the potential of ecosystems to 
provide more regulating services even compared to the control site in the study area. Therefore, 
plantations were able to enhance the potential of preserving the future of ecosystems and sustain 
the supply of ecosystem services. Fan et al. (2016a) concluded that Haloxylon sp. can restore 
arid ecosystems by increasing herbaceous plant cover and soil nutrients. Thinning of H. 
ammodendron plantations is proposed to reduce their excessive consumption of soil water (Zhu 
and Jia, 2011). Restoration of species composition and diversity also plays an important role for 
ecological success, and our results showed that H. ammodendron plantations are not successful 
in improving species composition and diversity. Hence, the plantation of H. ammodendron is 
relatively success ecologically. Gann et al. (2019) assessed principles for ecological restoration 
and indicated that ecosystem restoration is a means of improving ecosystems and human 
well-being in the social system. Improving the provisioning services that are directly related to 
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human well-being is most important for local people in the social system (Millennium Ecosystem 
Assessment, 2005). In this study, water yield and stocking rate at the 15- and 30-year-old 
plantation sites that have direct impacts on social welfare, respectively increased by 2.17 and 
5.43 times compared to the degraded site, but they failed to compete over the control site. 
Therefore, the plantation of H. ammodendron cannot be considered successful due to local 
conservation motivation in the social system. In general, monoculture of H. ammodendron will 
not be successful ecologically and socially. It should be suggested that the combined plantations 
of H. ammodendron with some other plant species that are compatible with H. ammodendron can 
further contribute to the success of restoration in arid ecosystems. 


5 Conclusions 


The method presented in this study helps managers to evaluate the ecological-social success of 
restoration in order to select restoration practices, especially plantation with non-native species in 
arid and semi-arid areas. Land cover changes resulted from plantations first impacted on 
ecosystem structures. In arid and semi-arid areas that are faced with environmental disturbances 
such as wind storms, the vertical structure of ecosystems is also an important parameter in 
addition to the horizontal structure. However, larger species are mostly preferred in plantations in 
arid and semi-arid areas due to their influences on improving the vertical structure of ecosystems. 
The results of this study showed that in arid and semi-arid areas, species diversity as an important 
attribute of ecosystem structures is more effective on improving ecosystem services than vertical 
structure, which should be considered in management plans in the future. In the process of 
restoring arid ecosystems, not only the vegetation structure and species interaction ecologically, 
but also ecosystem services socially should be considered as criteria for species selection. The 
plantation of non-native species that is structurally different from native species may change all 
the complex relationships between ecosystem structures and ecosystem services, so this should be 
considered in sustainable management. Compared to the degraded site, improving the economic 
values of ecosystem services at the plantation sites does not indicate the success of restoration. 
The measure of the success should be the degree to which a restoration practice meets the 
potential of the control site. Sustainable ecosystem management should not only consider 
ecosystem services related to local conservation to motivate local people to participate in the 
implementation of restoration projects and protect ecosystems, but also consider ecosystem 
services related to preserving the future to promote the ecological sustainability of ecosystems. 
The method presented in this study helps managers to evaluate the ecological-social success of 
restoration (especially plantations with non-native species in arid and semi-arid areas) to further 
select restoration practices. 
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Appendix 
The detailed information on the quantification of the seven ecosystem services is as follows. 


1 Soil stability 


To measure soil stability, we derived soil stability index (%) using the Landscape Function 
Analysis (LFA) methodology. For each patch identified along the transect, seven soil surface 
indicators were assessed: soil cover, litter cover, cryptogam cover, crust broken-ness, erosion type 
and severity, deposited materials and surface resistance to disturb and slake (Tongway and 
Hindley, 2003). The contingent valuation method (CVM) was used to estimate the economic 
value of soil stability. 

Supporting services indirectly benefit people by supporting other services (Millennium 
Ecosystem Assessment, 2005). Therefore, the indirect benefits of soil stability (for services such 
as water and livestock production) were explained by the respondents connected with nature and 
living around the sampling sites. The field investigation was conducted on May and June of 2022, 
and respondents (180 individuals) were asked "Do you agree with the benefits demonstrated and 
how much are you willing to pay for enhancing these benefits (supplying of soil stability)?". 


2 Aesthetic value 


Two 100-m transects were laid out at each sampling site. Thirty 2 mx2 m quadrats were randomly 
established along each transect. In each quadrat, canopy cover (%), litter cover (%), and stone and 
gravel percentage (%) were estimated. The number of plant individuals was also recorded. 
Previous researches showed that people's perception of aesthetic value can be enhanced by the 
abundance of flowering plants, as well as their colorfulness and blooming seasons (Graves et al., 
2017a, b; Uchida et al., 2020). Using the following equation, the aesthetic value for each quadrat 
was estimated (Schirpke et al., 2017): 

AV =rD+rFLP+rFCA+rFPCC , (1) 
where AV is the aesthetic value; rD is the relative Simpson's diversity index; rFLP is the relative 
flowering period; rFCA is the relative amount of flower colors; and rF PCC is the relative canopy 
cover of flowering plants. The parameter of rD was estimated by dividing the Simpson's diversity 
index of each quadrat by the maximum Simpson's diversity index of all quadrats; rFLP was 
obtained by dividing the total flowering time of the species in each quadrat by the total flowering 
time of all existing species in all quadrats; FCA was estimated by dividing the maximum number 
of flower colors observed in each quadrat by the number of flower colors in all quadrats; and 
rFPCC was calculated by dividing the canopy cover of flowering plants in each quadrat by the 
maximum canopy cover of flowering plants in all quadrats (Schirpke et al., 2017). 

We determined the flowering period of plants according to literature (Karimi et al., 2008; 
Toopchi-Khosroshahi and Lotfalizadeh, 2011; Ariapour et al., 2015). We also checked the flower 
colors observed in the quadrats based on botanical literature (Rechinger, 1997). CVM was used to 
estimate the economic value of aesthetic value. A sample size of 180 individuals was obtained 
using Cochran method (Cochran, 1997). Respondents (180 individuals) were randomly selected 
from visitors and people who are connected with nature and lived around the sampling sites. A 
questionnaire method was applied to assess people's willingness to pay for aesthetic value. 
Respondents were asked "How much are you willing to pay for a day or more on this site if you 
enjoy these views" based on the method of Garcia-Llorente et al. (2012). 


3 Carbon sequestration 


The Integrated Valuation of Ecosystem Service and Tradeoff (InVEST) model was used to 
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measure carbon sequestration. In this model, aboveground biomass carbon (Mg/hm/’), 
underground biomass carbon (Mg/hm7), litter carbon (Mg/hm’) and soil carbon (Mg/hm?) can be 
estimated. Thirty 2 mx2 m quadrats were randomly established at the control and degraded sites. 
For the two plantation sites, thirty 10 mx10 m quadrats were randomly established. In each 
quadrat, the canopy cover (%), height (m) and number of individuals of plant species were 
determined. Double sampling method was used to estimate the aboveground biomass (Mg/hm?) of 
the species (Reid et al., 1990). Hence, the aboveground biomass was visually estimated in each 
quadrat. However, a quarter of the quadrats were considered as calibrated quadrats, in which the 
aboveground biomass was also measured by direct harvesting (clipping and weighing) after 
estimation. The clipped samples were transferred to the laboratory to determine the moisture 
content and carbon content. The root/stem ratio (the ratio of the underground biomass to the 
aboveground biomass) was applied to estimate the underground biomass of the species. For this 
purpose, five average-size individuals were selected for each plant species and their root samples 
were taken by digging a soil pit to the rooting depth, and then the samples were weighted. A total 
of 100 g of root samples were taken to determine the moisture content and carbon content. Ten 
soil samples were randomly taken from 0-30 cm soil depths at each sampling site. The organic 
matter of all samples (i.e., aboveground and belowground biomass, litter and soil organic matter) 
was determined by wet and dry combustion. Finally, 11.00 USD was considered as carbon price 
per ton to estimate the economic value of carbon sequestration (Dang et al., 2022). 


4 Water yield 


Water yield was quantified using the water yield model in the INVEST model. The water yield 
model has been applied in different regions of Europe (Bangash et al., 2013), America (Hamel et 
al., 2015), Africa (Belete et al., 2020) and Iran (Daneshi et al., 2021; Sharafatmandrad and 
Khosravi Mashizi, 2021; Balist et al., 2022; Darvishi and Yousefi, 2022) to quantify water yield. 

In this model, water yield was estimated by subtracting evapotranspiration from precipitation in 
the region (Tallis et al., 2011). Required data in this model are as follows: precipitation (mm), 
reference evapotranspiration (mm), plant water availability (cm?/cm3), soil depth (m), crop 
coefficient and maximum rooting depth (m). The Penman-Monteith equation was used to estimate 
reference evapotranspiration based on daily mean, maximum and minimum temperatures obtained 
from local weather stations. Plant water availability was determined using soil texture data (De 
Ridder and van Keulen, 1995). Soil depth was estimated using the methodology of Tsai et al. 
(2001), Tesfa et al. (2009) and Liu et al. (2013) based on topography characteristics determined 
by digital elevation model (DEM) and soil cover. Crop coefficient was estimated according to the 
method of Allen-Wardell et al. (1998) and leaf area index was determined using the gravimetric 
method (Jonckheere et al., 2004). The maximum rooting depth was calculated according to the 
studies of Canadell et al. (1996), Schenk and Jackson (2002) and Tallis et al. (2011). Water price 
per cubic meters was considered to be 0.60 USD (Mousavi et al., 2021). 


5 Soil protection 


The InVEST model was used to estimate soil protection. This model can be applied to reveal the 
impacts of land cover change on soil protection (Hamel et al., 2015), which has been tested in 
many countries and regions of the world, e.g., Ethiopia (Aneseyee et al., 2019), China (Zhou et al., 
2019), Tana watershed, Kenya (Ayana et al., 2017), USA (Hamel et al., 2015) and Iran (Khosravi 
Mashizi et al., 2019). In this model, soil loss can be estimated using the universal soil loss equation 
(USLE) (Renard et al., 1997). Soil protection (S; Mg/(hm?-a)) was estimated by soil loss from 
non-vegetated areas (Slmax; Mg/(hm?-a)) minus the actual soil loss (SL;; Mg/(hm?.a)) as follows: 
S = Slmax — SLj 5 (2) 
Slmaxx=RXKXLXSXP, (3) 
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SL; =RXKXLXSXCxP, (4) 
where R (MJ-mm/(hm7-h-a)) is the total rainfall erosivity in one year; K (t-hm?/(MJ-mm)) is the 
soil erodiblity factor; Z and S are the slope length and steepness, respectively; P is the erosion 
control practice factor; and C is the vegetation management factor. R was estimated using the 
method of Wischmeier and Smith (1978). K was calculated using the method of Romkens et al. 
(1986) based on soil texture. Soil samples were randomly taken from 0-25 cm soil depths at each 
sampling site to determine soil texture. The hydrometer method was used to determine the particle 
size distribution (sand, silt and clay contents) of each soil sample (Soil Survey Staff, 1994). L and 
S were estimated using the DEM data based on the method of McCool et al. (1987). P was 
estimated using the Wener method (Lufafa et al., 2003), and C was determined based on the study 
of Wischmeier and Smith (1978). Soil price per ton was considered to be 0.82 USD (Rastgar et 
al., 2016). 


6 Habitat provision 


To estimate the diversity of subsoil macrofauna, we established ten 25 cmx25 cmx25 cm plots 
under the dominant shrubs at each sampling site, and soil samples were taken (Palacios-Vargas et 
al., 2007; Sandhu et al., 2010; Ghaley et al., 2014). Macrofauna of each soil sample were 
separated by hand and killed by immersion in ethyl alcohol (70%). Macrofauna were identified at 
the family level. The diversity of understory plant species was also measured in 30 quadrats (50 
cmx50 cm) established under the dominant shrubs at each sampling site (Bueno et al., 2013). 
Simpson's diversity index was used to calculate the diversity of understory plant species and 
macrofauna. The average value of Simpson's diversity index for the diversity of understory plant 
species and macrofauna was considered as the potential of habitat provision. CVM was used to 
estimate the economic value of habitat provision. Habitat provision as a supporting service 
indirectly provides benefits for people. Like soil stability service, the benefits of habitat provision 
were explained and respondents (180 individuals) were asked "Do you agree with the benefits 
demonstrated and how much are you willing to pay for enhancing the benefits (supplying of 
habitat provision)?". 


7 Stocking rate 


Forage supply (stocking rate) is the main provisioning service of rangelands (Oñatibia et al., 2015). 
It is one of the main services of rangeland ecosystems, and is drastically depends on management 
plans (Maes et al., 2012). Metabolizable energy in forage and metabolic energy required by 
livestock were used to estimate the stocking rate at each sampling site (Coupland, 1992; Robles 
and Passera, 1995). Metabolizable energy in forage was determined using forage quantity and 
quality. Thirty 2 mx2 m quadrats were randomly placed at the control and degraded sites. For each 
of the two plantation sites, thirty 10 mx10 m quadrats were randomly established. The amount of 
forage produced in each quadrat was measured by double sampling method. Metabolizable energy 
of forage was also estimated by determining nitrogen content (%) and acid detergent fiber 
(ADF; %) of forage in each quadrat (Minson, 1987). Nitrogen content of samples was determined 
by Kjeldahl method (Saez-Plaza et al., 2013), and ADF was measured using the method of Van 
Soest et al. (1991). Daily required metabolic energy by livestock energy requirement was also 
estimated using the Ministry of Agriculture, Fisheries & Food (MAFF) method (MAFF, 1975). 
Goat (average weight of 30-35 kg) is the dominant livestock in the study area. The economic value 
of stocking rate for each sampling site was estimated using the market price of meat. 
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